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1 Motivation 



Irradiation effects in graphite were one major research area in the past, par- 
tially due to its application as a moderator in thermal nuclear reactors. Graphite 
is still a material of choice for nuclear applications due to its low cross-section 
for neutron absorption. The influence of irradiation damage produced by dif- 



ferent kinds of ions on several properties of graphite was reviewed by iKellv 



(119811 ) in chapter 7 of his b ook. Recently i rradiation effects in carbon nanos- 



tructures were reviewed by Banhart (1999). The effect of neutron irradiation 
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on the magnetic properties of graphite has been studied in the past and shows 
the expected results, i.e. the introduction of lattice defects by irradiation pro- 
duces a decr ease in the diamagnetism and an increase in the spin density 
f)KelM 1 19811) . We are not aware of any study made in the past on the effects 
of proton irradiation on the magnetic properties of graphite or carbon-based 
structures. 



When we started this research work on irradiation effects on the magnetism 
of graphite we had two reasons to choose protons as energetic particle. The 
first one is related to the analysis method called PIXE (Particle Induced X-ray 
Emission) that uses protons to get a map for all relevant impurity elements 
within a sample depth of 30 /im for a proton energy of ~ 2 MeV in carbon (see 
Sect. 2.2 for details). A systematic and full characterization of the magnetic 
impurity content in each of the samples, and after each treatment or handling 
(it makes no sense to start with a highly pure sample and then cut it with a 
steel knife afterwards), is of primary importance and absolutely necessary. 



The second reason was based on early reports on room-temperature ferromag 



netic behavior in some carbon-based structures (see references in iMakarova 
From those early works our attention was focused to the magnetic 
properties found in amorphous-like carbon prepared from different hydrogen- 
rich starting materials where an increase of the saturation ma gnetization with 



the hydroge n concentration in the starting material was found ([Murata. Ushiiima. Ueda. and KawagucI 



1991. 11992T ) . The origin for the magnetic ordering has been related to the mix- 



ture of carbon atoms with sp 2 and sp 3 bonds, which was predicted to reac h 
a magnetization higher than in pure Fe (|Ovchinnikov and Shamovskvl . 11991). 
Hydrogen, on the other hand, was a ssumed to have a role only in the formation 
of the amorphous carbon structure ([Murata. Ushiiima. Ueda. and Kawaguchi , 



1992). New theoretical predictions, however, show that hydrogenated graphite 



can display spontaneous magnetizat ion coming from different number s of mono- 
and dihydrogenated carbon atoms ( Kusakabe and Maruvama . 20031 ) . Sponta- 
neous magnetization may also appear in the case of monoh ydrogenated zigzag 
edges ( Fuiita. Wakabavashi. Nakada. and Kusakabe) . 1996| ) if the distance be- 
tween them is large or if they are not in parallel. 



The advantage of proton irradiation is twofold: it enables us to make an impu- 
rity analysis simultaneously to the implantation of hydrogen. In this chapter 
we will review the main effects obtained after proton irradiation in differ- 
ent carbon-based structures. This chapter is organized as follows. In the next 
section we provide the main characteristics of our irradiation facility. In sec- 
tion 2.2 we show an example of element analysis obtained in one of the graphite 
samples used for the irradiation studies. The irradiation effects are reviewed 
in section 3. This section is divided in two main subsections that describe the 
effects in oriented graphite and carbon-based thin films. In section 4 we discuss 
some of the effects observed after annealing the sample at high temperatures 
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in vacuum or after leaving it at room temperature for a long period of time. 

In this chapter we concentrate ourselves mainly on the effects produced by 
proton irradiation in highly oriented pyrolitic graphite (HOPG). Effects of 
irradiation with alpha particles will be discussed shortly in section 3.1.3. Iron 
and ffuor irradiation effects on HOPG and on diamond are currently under 
way and the results will be published elsewhere in the future. 



2 Irradiation Characteristics 



All the irradiations presented in this chapter were done with the high-energy 
nanoprobe LIPSION of the University of Leipzig. The accelerator is a single 
ended 3 MV SINGLETRON™ with an RF-ion source for protons and alpha 
particles. The focusing system can deliver proton beams with diameters as low 
as 40 nm at very low currents of the order of 0.1 fA. For the present irradiations 
we worked at 2.25 MeV and currents up to 500 pA in the case of proton beam 
of diameter of 1 — 2 fim (microbeam). We have also irradiated samples with 
a broad beam of 0.8 mm diameter, energy of 2 MeV and currents between 
~ 50 and ~ 150 nA. The X-ray detector is an Ortec HPGe IGLET-X and 
subtends 187 msrad solid angle. The nume rical simulations di scussed in this 
section are based on the code SRIM2003 (|Ziederl . Il977-1985h . The samples 
were mostly attached to a Si substrate using a small amount of Varnish or 
a mixture of Varnish with ultra-pure graphite powder. In one case we fixed 
the sample on a messing plate with a hole to enable the free penetration of 
the energetic particles (H + and He + ). The magnetic moment of the substrates 
was always measured before attaching the sample. All the irradiations were 
done at nominally room temperature; the temperature of the sample was not 
controlled during irradiation. 



2. 1 The interaction of Me V-protons with graphite 



Figure 1(a) shows a typical profile for the number of defects created in a carbon 
matrix per 2.25 MeV-proton and per 1 /im interval. The penetration range of 
about 46 /iin is clearly visible. For samples thicker than this range one expects 
that the majority of the protons will be stopped, whereas for samples substan- 
tially thinner essentially all protons traverse the sample. This means that for 
20 //m or thinner samples the implantation of protons should be undetectable. 

Strikingly, en hancement of ferromagnetism has also been observed in disor- 

dered carbon ( Hohne. Esquinazi. Han. Spemann. Setzer. SchaufuB. Riede. Butz. Streubel. and Hessel . 
12004 ) and fullerene films deposited on Si substrates after proton irradiation 
(see section 3.2). Further details on the different processes that influence the 
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Fig. 1. SRIM2003 Monte Carlo simulations (jZiegleTJ . Il977-1985h for 2.25 MeV pro- 
tons on HOPG. (a) The number of defects per ion and 1 /im depth interval. The 
penetration range of 46 /jm is clearly visible, (b) Sketch of the modified area in 
graphite due to 2.25 MeV proton bombardment. The lateral straggling has a full 
width at half maximum (FWHM) of ~ 2.5 fim for an initial beam diameter of 1 /im. 
The numbers are obtained as suming a displacement energy of 35 eV for Frenkel 
pairs in HOPG. Adapted from lButz. Soemann. Han. Hohne. Setzer. and Esauinazil 



penetration of protons in gra phite as well as the restrictions of the SRIM2003 
simula tions were discussed bv lButz. Spemann. Han. Hohne. Setzer. and Esquinazi 



The question how many protons remain in thinner samples cannot 
be clearly answered yet. Experiments with thinner and unsupported samples 
(to rule out part of the backscattering processes) are necessary and will be 
performed in the future. Under the assumption that there is no annealing 
of defects during irradiation and that the damaged area is of the same size 
as the irradiated one, the numerical simulations indicate that for fluences 
0.001 ... 75 nC//zm 2 we get in the near surface region between 4.7 x 10~ 6 and 
0.35 displacements per carbon atom, i.e. complete amorphization for the high- 
est fluence, using a displacement energy of 35 eV for Frenkel pairs in HOPG, 
in agreement with re cently published studies of the damage casca des by irra- 
diation on graphite ((Abe. Naramoto. Iwase. and Kinoshital . 119971 ) . Changing 
this number to 28 eV yields defect densities which are about 30% higher. To- 
wards the end of the track, the lateral beam straggling becomes important and 
the damaged area can be several microns wide, see Fig. 1(b). For a fluence of 
75 nC/nm 2 we have ~ 5 x 10 11 protons//zm 2 , i.e. the regions where defects 
are created by each individual proton overlap heavily. 



Note that even for the highest fluence the relative concentration of hydrogen 
(H/C, see Fig. 1) in the first micrometer from the surface remains rather low, 
of the order of 10 -4 or lower. However, most of what happens after the defect 
formation is still unclear. A dangling bond could attract a hydrogen atom - not 
necessarily from the proton implantation but already present as impurity in 
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the sample. Interestingly enough, there is little information on residual hydro- 
gen from pyrolysis of hydrocarbons, possibly because detection methods for 
hydrogen contents of the order of 100 ppm or below were not readily available 
and such quantities were of little relevance for applications as moderators in re- 
actors. Hydrogen atoms in the van der Waals gap or within the graphite layers 
should be highly mobile, contrary to hydrogen atoms trapped at defects. In- 
vestig ations with deuterium done bv ISiegele. Roth. Scherzer. and Pennvcook 
(1993|) indicate that deuterium does not readily diffuse out of HOPG but is 
rather chemically bound up to D/C ratios of about 0.45, probably at lattice 
defects. The maximum retention of deuterium depends on the temperature 
and implantation energy. A broad range o f binding energ ies for hydrogen in 
graphite up to about 4 eV was reported bv lAtsumil ( 20021 ). i.e. even tempera- 
tures as high as 3000-3500 K may be not sufficient to eliminate all incorporated 
hydrogen. On the other hand, the effective activation energies for hydrogen 
diffusion in our system are not necessarily the same as, for example, those 
obtained by experimental methods - usually at high temperatures - to study 
kinetics of diffusion of hydrogen in graphite. Irradiation effects of M e V pro - 
tons on diamond-like films were studied by IWang. Wang, and Chenl ([1992). 
According to the authors, below a fluence of the order of 1 nC//zm 2 the hy- 
drogen atoms produced by ion irradiation could be recaptured by dangling 
bonds. For higher ffuences a release of hydrogen is expected. 



Irradiation at a fixed energy has a clear disadvantage. If the magnetic ordering 
is triggered at a specific density of protons and defects, then, it is clear that 
with fixed proton and defect distributions, such as the one shown in Fig. 1, 
we would have a rather narrow window to get a maximum effect. Therefore, 
it should not be a surprise if for some irradiation energies and fluences one 
measures negligible effects or even a reduction of the magnetic order present in 
the sample before or after some irradiation steps. Carbon magnetism as well 
as irradiation effects on it belong to a new field in magnetism. For graphite 
bulk and thin film samples there are still many questions to be clarified in the 
future concerning the hydrogen implantation by irradiation and its effects on 
the magnetism. 



The defect formation process by high energy protons is a non-equilibrium 
athermal process and it appears rather unlikely that ordered arrays of defects 
are formed by migration of interstitial carbon atoms or vacanc ies, maybe with 
the exception of the interstitial across the gallery. According to Banhart ( 19991 ) 
and from electron irradiation studies, the essential types of radiation damage 
up to intermediate temperatures are the rupture of basal planes (due to shift 
of the C-atoms out of the plane) and the aggregation of interstitials into 
small dislocation loops between the graphene layers. The migration energy of 
the interstitial depends whether it is bounded. Di-interstitials were proposed 
to explain the irra diation- i nduce d amorphization of graphite with a migration 
energy of 0.86 eV ( Niwasel . 1995| ). The interstitial loops are stable up to rather 



5 



o 


E 


CO 


Oc 


>, 


oo 

LO 


'</> 








B 






v ~ 0.1 : 




E 


CD 


o 


Ram 


1360 



0.01 



0.1 




Fluence (nC/(im ) 



c 
=i 
>* 

1— 
CO 
1— 

-2- 





\ (a) 

I 




1 (b) 





290 288 286 284 282 
Binding Energy (eV) 



Fig. 2. Left figure: Raman intensity ratio as a function of the irradiation fluence 
for spots created on a HOPG sample (see section 3.1.2). The line represents the fit 
curve O.lx 0,5 . Right figure: Decomposition of the XPS C\ s core-level peaks of HOPG 
surface after (a) and before proton irradiation (b). The components after peak fit 
are: at 284.4 eV: "C-sp 2 "; at 285.0 eV: defect peak "C-sp 3 "; other components: 
different C-0 bonds. The thick and thin solid lines denote the experimental curves 
and the components, respectively. The measurement was done after irradiation of 
the sample with 1.2 mC and a fluence of 0.6 nC/mm 2 (ste p 4 of t he sample reported 
by Esauinazi. Soemann. Hohne. Setzer. Han, and Butzl d2008l)). Adapted from 



Hohne. Esauinazi. Han. Soemann. Setzer. Schaufufi. Riede. Butz. Streubel. and Hessel 
2C 



1993). Certainly, the layer 



high temperatures, probably to 1000° C (jBanhart . 
structure should remain essentially intact, at least for our low fluences and at 
the sample surface. 

Irradiation can also produce a transition from sp 2 to sp 3 bonding leading 
to cross- l inks b etween the graphene layers and the formation of sp 3 clusters 
( Tanabel .ll99rl These clusters appear to be stable and do not anneal at high 
temperatures. Our own Raman and X-ray Photoemission Spectroscopy (XPS) 
measurements indicate the following. For all unirradiated areas of the HOPG 
samples the Raman spectra show only one pronounced peak at 1580 cm -1 
(the E 2 g2 mode) as expected for graphite without any disorder. For the spots 
produced with the proton microbeam the Raman spectra show additionally 
the disorder mode D at 1360 cm -1 . With increasing fluence the E 2 g 2 and D 
modes become broader and the ratio of their intensities I(D)/I(E 2 g2), as a 
measure of the degree of disorder, increases (see Fig. 2). 



The XPS results show a clear difference between an irradiated and an unir- 
radiated sample, see Fig. 2. The Ci s core-level peaks of HOPG samples were 
recorded using a pass energy of 10 eV and were fitted by six components with 
binding energies of 284.4 eV (main peak of HOPG, C-sp 2 ), 285.0 eV (defect 
peak "C-sp 3 "), 285.8 eV, 286.3 eV, 289.0 eV (C-0 components) and 290.9 eV 
(shake-up peak). The main result is that the defect peak component increases 
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after proton irradiation. A similar behavi our was found for HOPG under 



plasm a low-energy argon-ion bombardment ([Rousseau. Estrade-Szwarckopf. Thomann. and Brault , 

Whether this defect peak corresponds to a pure C-sp 3 state or whether 
some of the "defect" carbon atoms are bounded to hydr ogen is unclear and it is 
a matter of current research ((Estrade-Szwarckopa . 120041 ) . The knowledge of the 
C-H sp 2 and sp 3 bonds is for the magnetic properties of carbon structures of 
importance. We remark that numerical simulations indicate a 100% polarized 
7r-band, i.e. a ferromagnetic order stable at room temperature, for a graphene 
layer with a mixture of sp 2 and sp 3 b ondings (mono- and di-hydrogena ted) at 
the zigzag edges of a graphene layer ( Kusakabe and Maruvama , l2003i l . 



Other aspect of the irradiation effects on graphite that may be of importance 
to trigger magnetic ordering is the formation of pentagons and hep tagons in 



the b asal planes, which could cause a bending of the graphene layers ([Banhart 



1999). The influence on the electronic band str ucture of graphite of such topo- 



logica l defects has been studied theoretically by Gonzalez. Guinea, and Vozmedianol 
( 20011 1 . According to these authors, these defects may trigger ferromagnetism 
or even superconductivity. 



For thick targets and sufficiently high ffuences, the target surface actually 

swells, which can be easily measured by an atomic force microscope ( AFM) , see 

Fig. 8(a), or is even visible under an optical microscope (Spemann. Han. Esquinazi. Hohne. and Butzl . 



2004 : Esquinazi. Han. Hohne. Spemann. Setzer. and Butz . 120051 ). The swelling 
in the c— direction occurs together with the contraction in the graphene layer; 
the new formed interstitial planes push the existing planes apart leading to 

a protuberance at the sample surface. This i r radia ti on effect was studied by 

Koike and Pedraza ( 1994 ): Muto and Tanabel ( 1997 ) : Takeuchi. Muto. Tanabe. Arai. and Kurovanag 



(119971 ). This effect was not observed for thin enough targets for the usual pro- 
ton irradiation ffuences used in this work. 



Another important aspect regarding the irradiation effects is the heat load due 
to the beam. For a spot of 2 /mi x2 /im and a current of 500 pA this amounts 
to ~ 300 W/mm 2 for thick samples. Most of the energy will be deposited near 
the end of the ion track. Depending on the lateral dimensions and the thickness 
of the sample as well as the target holder, heat is transported away from the 
beam spot. The details are difficult to calculate, but the measurement and 
control of the temperature is of importance. For thin samples, this is relatively 
unimportant. We have evidence that at higher currents the magnetic response 
in MFM-measurements is strongest at the rim of the spot suggesting annealing 
effects in the center of the spot, as Fig. 3 indicates. For lower currents the 
magnetic image is homogeneously distributed across the spot (see section 3.1), 
in contrast to the magnetic "ring" shown in Fig. 3. 
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Fig. 3. (A) Topography and (B) MFM images obtained from a spot produced with 
a dose of 6.8 nC/^m 2 and a current / = 0.51 nA, which means a current density 
larger than 5 mA/cm 2 . The scan area is 15 /imxl5/im. The bottom figures show the 
line scans done at the positions of the black lines of the upper pictures. The MFM 
images were taken at a distance of 50 nm using standard m agnetic tips. Taken from 
Esauinazi. Han. Hohne. Spemann. Setzer. and Butzl (|2005h 

2.2 Impurity measurements 



With the proton micro-beam a total charge of 0.5 fiC suffices to obtain a 
minimum detection limit for Fe impurities below 1 yUg/g. With our broad- 
beam PIXE facility we work with 2 MeV protons and currents of 150 nA 
and a beam diameter of 0.8 mm. Here, about 5 /itC are required to reach 
a minimum detection limit of 1 fig/g for Fe. The fluences are in the range 
of 1 . . . 50 /iC/mm 2 for the micro-beam and the broad-beam, respectively. 
These fluences are several orders of magnitude lower than the fluences ac- 
tually applied for the induction of magnetic ordering in graphite, as will be 
discussed below. The advantage of the micro-beam compared to the broad- 
beam is that a distribution map for all relevant impurity elements can be 
obtained contrary to the integral value for the broad-beam method, a rather 
important issue considering the grossly inhomogeneous distribution of impu- 
rities like Fe which we have encountered in some samples (see for exa mple 
Fig. 1 in ISpemann. Han. Hohne. Makarova. Esquinazi. and Butzl fj2003la A 



typical broad beam PIXE spectrum for a HOPG sample is shown in Fig. 4. 
It shows the presence of a number of impurities, the Fe content being (0.45 ± 
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Fig. 4. Typical broad beam PIXE spectrum from a HOPG sample. The main im- 
purities are Ti, V, and Fe. The minimum detection limit for other elements heavier 
than Si is ~ 0.3 /xg/g. Note that 2.1 /ig/g iron in the carbon matrix means a con- 
centration of 0.45 ppm Fe. The background curve is due to the carbon matrix. 



0.04) ppm. From the study of ferromagnetic signals in different HOPG samples 

( Esauinazi. Setzer. Holme. Semmelhack. Kopelevich. Spemann. Butz. Kohlstrunk. and Loschd . 



2j) we may suspect that HOPG could contain a non negligible amount of hy- 
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drogen, some of them may be related to the origin of the magnetic signals. The 
amount of hydrogen before or after irradiation within the fruence used in this 
work cannot be determ ined by PIXE. Taking into account recently done the- 



oretic al work done by iLehtinen. Foster. Ma. Krasheninnikov. and Nieminen 



()2004f ) that estimates a magnetic moment of ~ 1 hb for a carbon vacancy, 
2.3 [is in the surrounding of an hydrogen bonded to a carbon atom at the 
position of a carbon vacancy, or 1.2 n B for a vacancy that is saturated by two 
hydrogen atoms in a graphene layer, we may expect that a few ppm hydrogen 
trigger a non negligible magnetic signal in graphite. Therefore, hydrogen con- 
tent measurements should be carried ou t with a sensitivity in the ppm range, a 



rather difficult task. The reader can read lButz. Spemann. Han. Hohne. Setzer. and Esquinazi 



(2005) :or a short discussion on the methods available for this kind of mea- 



surement. 
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Fig. 5. Magnetic force microscope (MFM) and SQUID measurements of a HOPG 
sample irradiated with 4 broad spots (</> ~ 0.8 mm) of 135 //C total charge each. All 
images scan an area of 20 x 20/um 2 : (a) MFM image of the sample in the state just 
after irradiation, (b) after application of a field of + 1 kOe, (c) after application 
of a field of -1 kOe. All measurements were done in remanent state after removing 
the field, i.e. at no applied field, (d) Magnetic domain structure of the sample in 
the same state as in (a) but at a different position. All the MFM measurements 
were done with a standard magnetized MESP tip and at a distance of 50 nm from 
the sample surface, (e) The magnetic moment of the same sample measured with 
a SQUID at room temperature. The sample was fixed on a Si substrate for the 
irradiation. The data in the figure are the original data of the sample with Si 
substrate, without any substraction. The signal of the sample overwhelms that from 
the substrate. Measurements of the sample standing alone reproduces basically the 
curve presented here. 
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3 Irradiation Effects 



3. 1 On highly oriented pyrolytic graphite 
3.1.1 Broad proton irradiation 

All the irradiated HOPG samples discussed in this review were from Advanced 
Ceramics Co. (ZYA grade, 0.4° rocking curve width at half maximum) with 
a content of magnetic metallic impurities below 1 ppm. The largest concen- 
tration of non magnetic metallic ions was found for Ti (< 6 fig/g) and V 
(< 16 /ig/g). A typical surface area of the samples was 2x2 mm 2 and a 
thickness between 0.1 mm to 0.3 mm. The samples were glued with varnish 
(or a mixture of varnish and a high purity graphite powder to increase the 
thermal coupling) on a high-purity Si substrate and the magnetic moment 
of the whole ensemble as well as of the Si substrate alone was measured. In 
general and within experimental error the diamagnetic hysteresis loops for the 
used Si substrates are reversible; after subtracting the linear diamagnetic de- 
pendence there is no significant hysteresis left (for an example see Fig. 7(a)). 
Nevertheless, when the ferromagnetic moment of the sample is weak, it can 
be a hard task to obtain the true sample magnetic moment from the SQUID 
signal of the mixture (Si plus the HOPG sample) and several checks have to 
be done. 

In this section we will discuss results obtained after irradiation of a large area 
of the HOPG sample using the micro- and the broad-beam of protons. As we 
will realize below, there are several irradiation parameters that may have an 
important role in inducing the achieved magnetic signal. Namely: 

- (1) The total implanted charge. The total amount of protons or alpha par- 
ticles that were implanted in the sample or travelled through the sample. 

- (2) The input energy. Although with the LIPSION we have the possibility of 
changing this energy from 1 to 2.25 MeV, in order to minimize the number of 
variables, the studies presented here were done with fixed irradiation energies 
of 2.0 or 2.25 MeV. This means that we have a well defined defect and im- 
plantation profile inside the sample, see Fig. 1. If the highest magnetic signal 
is determined at a given defect and proton density, it is clear that successive 
irradiation at similar energies may produce contrary effects if this optimized 
region is destroyed by the next irradiation at similar conditions. 

- (3) The fluence, the irradiated charge per unit area. 

- (4) The ion current. Large currents might heat the sample and non system- 
atic effects are then possible. 

- (5) Micro- or macro-irradiation, i.e. a broad irradiation with the 0.8 mm 
beam or a large number of micrometer spots distributed in the sample. Ex- 
perience indicates that the magnetic signals are much larger when one uses 
a high density of micrometer spots rather than the implantation of similar 
amount of charge with a broad beam. This fact appears to be related with the 



11 



density and/or type of defects that the beam produces. 

- (6) Sample temperature. All the irradiations presented in this review were 
done at nominally room temperature. Future experiments should try to irra- 
diate at lower or higher temperatures to check for its influence. 

- (7) Finally, the initial state of the sample, namely the density and type of 
defects, which in part can determine its metallicity. Computer simulation re- 
sults of the effects of adsorbed hydrogen on the band structure of a graphene 
layer indicate t hat metallization caused by specific de fects can quench a spin 
polarized state ( Duplock. Scheffler. and Lindanl . 2004 ). 



The first SQUID measureme nts that indicate a magnetic ordering after proton 
i rradia tion were published bv lEsquinazi. Spemann. Hohne. Setzer. Han, and Butz 



(2003). n that work an increase in the hysteresis loop was observed after sev- 



eral irradiation steps. These irradiation steps contained successively irradia- 
tion of several thousands of micrometer small spots as well as four (or three) 
spots of 0.8 mm diameter each on the same sample. In what follows we shall 
discuss further work that has been done afterwards, where we have tried to 
characterize the effects produced by some of the irradiation parameters de- 
scribed above. We stress that they are the very first steps to get reproducible 
magnetic order in carbon, a task that turned to be full of difficulties, as the 
published results from literature indicate. 

A broad and homogeneous proton irradiation of usual fluences (~ 150 iiC 
total charge in an area of ~ 1 mm 2 ) per spot on oriented graphite produces 
a magnetic signal that in general can be well observed with a magnetic force 
microscope (MFM) on the irradiated surface. Figure 5 shows the MFM im- 
ages obtained on an irradiated area with a total charge of 540 fj,C distributed 
into 4 spots, 135 /iC each. After this broad irradiation there is no significant 
change in the topography worth to note. However, the periodic magnetic do- 
mains are well defined. In the remanent state just after irradiation, in some 
of the irradiated area of the sample we observed two domain structures, one 
normal to the other, see Fig. 5(a). The period of the small domain structure 
is ~ 0.8... 1.2 fim depending on the region (compare with (d) which was 
obtained in the same state but in other irradiated area), whereas the other 
domain structure has a period of ~ 10 //m and depends on the scan direction of 
the MFM tip. After application of a field of 1 kOe in the +z direction (perpen- 
dicular to the graphene layers) the "vertical" domain structure shows a period 
of ~ 1.6 . . . 2.0 /zm and the width of the domains increases, see Fig. 5(b). The 
other domain structure is not observed. After magnetizing the sample with a 
field in the other direction there is a slight change of the width of the bright 
relative to the dark regions, compare (c) with (b). In non- irradiated graphite 
areas we did not find any signature of domain structures within the resolution 
of the microscope. This irradiation triggered a relatively large hysteresis loop 
that could be very well measured with the SQUID without any background 
substr action, see Fig. 5(e). The temperature dependence of this signal was 
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studied removing the sample from the Si substrate. It shows a weak decrease 
(~ 5%) with temperature of the saturation magnetization and coercivity fields 
between 5 K and 300 K. 

To study the influence of some of the irradiation parameters on the magnetic 
response of HOPG samples, we have studied three HOPG samples from the 
same batch. One sample of mass m = 1.19 mg (labelled AS 171 103/2) was 
irradiated with the proton microbeam as follows. First irradiation consisted 
of 10 4 spots of diameter (ft ~ 1.8 ^m each with a charge of 0.85 nC (total 
charge 8.5 /iC), fluence = 0.34 nC//im 2 , irradiated area: 0.62 x 0.62 mm 2 and 
current I ~ 1.1 nA. Second irradiation produced the same amount of spots but 
with a diameter of <ft ~ 2 /im and 0.72 nC charge each (total charge 7.2 /zC), 
fluence = 0.24 nC//iin 2 and / ~ 0.9 nA. The third and last irradiation was 
identical to the first one. The room temperature SQUID results of the sample 
before and after substraction of the diamagnetic background are shown in 
Fig. 6. It is clearly seen that the irradiation increases the magnetic signal 
of the sample. The last irradiation, however, did not produce any significant 
change respect to the last one. The remanent magnetization changed from 
M r (B = 0) ~ 3.3 x 10~ 5 emu/g to 2.7 x 10~ 4 emu/g and 4.5 x 10~ 4 emu/g 
for the first, virgin state and second (or third) irradiation steps. 

To check the influence of the irradiation current under broad irradiation condi- 
tions we have irradiated two samples with the following characteristics. Irradia- 
tion with a current / = 54±1 nA: sample AS211003/1, mass m = 1.33 mg, first 
irradiation: 4 spots of <ft = 0.8 mm and 53.4 /zC each (total charge 214 /xC); sec- 
ond irradiation: 4 spots of (ft = 0.8 mm and 110 /iC each (total charge 440 /xC); 
third irradiation: 4 spots of (ft = 0.8 mm and 54.5 /iC each (total charge 
218 /iC). The SQUID results after substraction of the diamagnetic background 
are shown in Fig. 7(a). We see that the initial magnetic state of this sample 
(magnetic moment at saturation m ~ 10 -6 emu) decreased after the first irra- 
diation. For the subsequent irradiation steps the magnetic moment increased. 
The saturation magnetization showed a decrease from its initial value of M r = 
(1.95 ±0.15)10~ 4 emu/g to M r = (1.0 ±0.08)10~ 4 emu/g after the first irradi- 
ation. For the 2 nd and 3 rd irradiation we have M r = (1.56 ± 0.15)10 -4 emu/g 
and M r = (2.6 ± 0.15)10~ 4 emu/g. The coercive field remains in the range of 
150 Oe. 

On the sample AS181103 of mass m = 1.12 mg, the following irradiation steps 
at a proton current of / = 150 nA were performed: (1) 2 spots of 0.8 mm diam- 
eter with 159 /iC each (total charge 318 /iC); (2) 2 spots of 0.8 mm diameter 
with 160 /iC each (total charge 320 /iC); (3) 2 spots of 0.8 mm diameter with 
154 /iC each (total charge 308 /iC). Figure 7(b) shows the SQUID results for 
this sample. Under those irradiation conditions the sample did not show any 
significant increase in the magnetic signal within experimental error. Com- 
paring this result with that of (a) we would conclude that the proton current 
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Fig. 6. (a) Magnetic moment vs. applied field of the HOPG sample AS171103/2 
glued on a Si substrate before and after irradiation (for details see text), (b) The 
same data as (a) but after substraction of the diamagnetic background. 
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Fig. 7. (a) Magnetic moment (minus the diamagnetic background) vs. field of a 
typical Si substrate (o) where all HOPG samples were attached before and after 
all irradiation steps. The other symbols correspond to sample AS2 11003/1 after 
different proton irradiations (see text for details) with proton current I ~ 55 nA. 
(b) The same for sample AS181103 after different proton irradiations with a current 
/ ~ 150 nA (for details see text). 
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plays an important role in inducing the magnetic ordering in HOPG samples. 



3.1.2 Magnetic spots and lines of micrometer size 



As mentioned in section 2, the LIPSION accelerator has the possibility to 
irradiate the samples with a proton micro- or nanobeam. We have used the 
microbeam to produce magnetic spots on oriented graphite surfaces. As for 
the broad beam irradiation, the microbeam was directed onto the HOPG sur- 
face parallel to the c— axis of the sample without beam scanning (excepting 
line scans, see below) leading to the formation of micron-sized spots with 
enhanced defect density, as measured by micro- Raman, see Fig. 2. In gen- 
eral, two spots separated by a distance of 20 /mi were irradiated with the 
same ion fluence and several ion fluences were used. For large enough fluences 

the swelling at the spots can be directly observed with an optical microscope 

(Soe mann. Han. Esquinazi. Hohne. and Butz, 2004; Es quinazi. Han. Hohne. Spemann. Setzer. and Bu 
12005b . The height of this swelling depends on the irradiated fluence and on the 



mass of the ions. The dependence of the maximum swelling height, measured 
by AFM, with the fluence for proton irradiation is shown in Fig. 8(a). With a 
MFM one can measure the magnetic signal on the spots. The maximum am- 
plitude of the signal (maximum phase shift) as a function of the fluence for two 
different proton currents is presented in Fig. 8(b). The units of the magnetic 

signal from the MFM are "degrees" ; a relat ion of this phase shift of the tip vi- 

bratio n to the force gradient can be found inlLohau. Kirsch. Carl. Dumpich. and Wassermann 



netic moment calcu lated from measurements of the phase shift in carbon sam- 

pies can b e found in Han. Soemann. Hohne. Setzer. Makarova. Esquinazi. and Butz 
f)2003bl lc]): Irian and Esauinazil (|2004h ~ 



We note that, whereas the swelling height increases systematically with flu- 
ence, the maximum phase shift of the magnetic signal does not, see Fig. 8(b). 
This result indicates that: (1) there is a negligible influence of the topography 
onto the magnetic signal and (2) the magnetic signal depends on the implanted 
charge, proton current and probably also on the carbon structure. The differ- 
ence between the behaviors observed for two different proton currents can be 
understood if we take into account that the higher the current the higher 
would be the internal temperature of the sample at the spot. Therefore, an- 
nealing effects can be the origin for this difference. It is worth to note that 
the signals could be reproduced after leaving the samples in air and at room 
temperature for two months, see Fig. 8(b). This is not the case for much larger 
times, see section 4. We note also that at the higher fluence range, the struc- 
ture of the spots is very probably that of highly disordered or even amorphous 
carbon. For both proton currents the maximum magnetic signal decreases at 
the highest fluences, see Fig. 8(b). 
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Fig. 8. (a) Maximum height of swelling as a function of the fluence, measured 
at micrometer small spots with an AFM. The three symbols indicate the height 
measured before and after the application of a magnetic field in z (up) and —z 
(down) direction. The measurements were done without applied field. The line is 
the function 27. 5x 0,61 nm with x being the fluence in nC///m 2 . (b) Maximum phase 
shift measured at micrometer small spots on two different HOPG samples, before 
application of a magnetic field, as a function of the irradiation fluence. The spot 
areas were 1 ^m 2 (B,*) and 4 /im 2 (A); the corresponding proton currents were 
171 pA and 855 pA, respectively. The (*) symbols correspond to the same sample 
measured again two months later. The lines are only a guide to the eye. 



Figure 9 shows an example of a magnetic spot as measured with the MFM be- 
fore and after application of a magnetic field in two directions. The irradiated 
region can be clearly recognized as the "white" spot of the magnetic images in 
(a). A line scan through these images indicates clear changes of the phase shift 
after application of a magnetic field, see top figure in Fig. 9(b). Usually the 
magnetic images are obtained with a tip to sample distance of 50 nm. If we 
increase this distance the phase shift amplitude decreases. This dependence 
can be used to estimate, although with a relatively large error, the order of 
magnitude of the magnetization at the spot surface. From the data shown in 
the bottom figure o f Fig. 9(b) one is ab l e to e stimate a magnetization of the 



order of 400 emu/g (H an and Esquinazil . 120 04). a very large magnetization of 



the order of soft ferromagnetic metals. 

As a first attempt to write magnetically on a graphite surface we have started 
with a simple cross. Figure 10(a-c) shows the magnetic images for two crosses 
made at different irradiation fluences and measured with different tips and the 
topography for one of them (d). The magnetic image changes after applying a 
magnetic field. The results after magnetizing the cross show a state with one or 
three magnetic domains. Micromagnet ic simulations with the appropriate pa- 



rame ters provide similar MFM images fS chmdler. Spemann. Ziese. Diaconu. Schmidt. Esquinazi. and I 



l2005h and they may be used to obtain material parameters (magnetic anisotropy, 
for example) that are not yet possible to measure directly. 
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Fig. 9. (a)Magnetic force gradient images (20 x 20 /xm 2 ) of a spot and its sur- 
roundings irradiated with 0.115 nC/um 2 (proton current I = 171 pA). The images 
were taken, from top to bottom, before field application, after applying a field of 
~ 1 kOe in the +z direction parallel to the c-axis, and in the —z direction. The 
tip-to-sample distance was 50 nm. (b) Top: The corresponding phase shift obtained 
at the line scans (white straight lines in figures (a)). The spot is located between 
~ 5 fim and ~ 10 am. The bottom, upper and middle lines in this region corre- 
spond to measuremen ts before and after application of a field in — z and +z d i- 
rection. Adapted from lHan. Spemann. Esquinazi. Hohne. Riede. and Butz 1 (2003a). 
Bottom: Scanning height dependence of the maximum phase shift at proton irradi- 
ated spots with fluences of 0.2 nC///m 2 and 2. 2 nC/ffm 2 . Solid lines are fits with 
the point probe approximations. Adapted from Han and Esquinazil (j2C 



3.1.3 Irradiation effects with alpha particles 



With the LIPSION accelerator we have also the possibility to irradiate the 
samples with a microbeam of 1.5 MeV alpha (He + ) particles. The SRIM2003 
simulation indicates that the penetration depth of these particles in graphite 
should be 4 /im. For a fluence of 0.1 nC//im 2 the defect density at the sur- 
face should be 1.4% and the region at the end of ion range should be rather 
amorphous. As for protons, we have produced pairs of spots at different flu- 
ences and studied the magnetic signal at the spot position with MFM. Figure 
11 shows the topography and magnetic images and the line scans of a spot 
produced with a fluence of 0.1 nC///m 2 . Whereas the swelling produced by 
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Fig. 10. (a) Magnetic image (top, 8 x 16 /Lim 2 ) and the corresponding phase shift 
(bottom) obtained at the line scan (white lines in all figures) for a cross produced 
with a fiuence of 0.5 nC/^m 2 and a proton current of / = 120 pA. The measure- 
ment was done before applying any magnetic field and with a LM (low moment), 
well magnetized tip. (b) Similar to (a) but for a cross produced with a fiuence of 
2 nC//xm 2 and after application of a field of 1 kOe normal to the surface. The 
magnetic image has a size of 18 x 18 /mi 2 and was taken with a standard MESP 
tip. (c) The same cross as in (b) but after application of a field of 2 kOe in the 
same direction as in (b). The magnetic image has a size of 10 x 20 /mi 2 , (d) Three 
dimensional image of the topography of the cross in (b,c) and the corresponding 
line scan (bottom). 
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Fig. 11. Micrometer spot produced after irradiation of alpha particles of 1.5 MeV on 
a HOPG surface with a fluence of 0.1 nC///m 2 . The upper pictures show the topog- 
raphy (left) and the magnetic images before application of a magnetic field (B=0), 
and after the application of a field in +z (B-up) and in — z direction (B-down). The 
bottom pictures show the line scans obtained at the position of the white lines in 
the upper images. 

the He + irradiation is remarkable, we recognize that the maximum phase shift 
measured by the MFM at the spot and its surroundings, see Fig. 11, is much 
smaller than those obtained at the spots produced with protons, see Fig. 9. 
Nevertheless, the fact that we can measure small but not zero signals, indi- 
cates that ion-, other than proton-, irradiation in the carbon structure may 
also trigger magnetic ordering. As we mentioned above, we do not know yet 
whether the hydrogen already present in the sample before irradiation or in 
the irradiation environment (rest gases in the chamber, for example) plays any 
role and therefore we should take these results with caution as a proof for a 
hydrogen-independent magnetic ordering. Future work should try to measure 
HOPG samples irradiated with alpha particles at different fluences and with 
the SQUID. 

3. 2 On thin films 

Taking into account the irradiation profile of protons in carbon obtained from 
Monte Carlo simulations, see Fig. 1, one would expect no or negligible effects 
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Fig. 12. (a) Sketch of the disordered carbon film CH577a produced by pulsed 
laser deposition with the irradiated (8 spots with 0.8 mm diameter and 150 /j,C 
total charge each) and non-irradiated areas, (b) Topography (AFM) and mag- 
netic (MFM) images in both areas, (c) Magnetic moment of the film (after 
substraction of the linear diamagnetic background) as a function of magnetic field 
applied parallel to the main area. The symbols (□,*) indicate: before irradiation, 
measurements done at T = 5 K and 300 K, and after proton irradiation, (•, *) 
at 5 K and 300 K. The inset shows an enhanced part around zero field. Adapted from 



Hol me. Esauinazi. Han. Spemann. Setzer. Schaufufi. Riede. Butz. Streubel. and Hesse] 
2C 



on the magnetic properties of samples, which thickness is much less than 
~ 40 fxm, because most of the protons would go through the sample leaving a 
small density of defects. We have checked this assumption irradiating carbon 
and fullerene films with thickness below 1 fxm. Figure 12 shows the MFM and 
SQUID results of a carbon film produced by pulsed laser deposition (PLD) in 
a hydrogen (H 2 ) atmosphere and deposited on a Si substrate. On the film 8 
spots of 0.8 mm diameter were irradiated (see sketch in the figure). The MFM 
signal shows magnetic domains in the irradiated region only, see Fig. 12(b). 
The SQUID shows a clear increase in the magnetic moment after irradiation, 
see Fig. 12(c). Because the results are taken on an disordered sample, they 
also reveal that a graphite ordered structure, at least in the mesoscopic level, 
is not necessary to trigger magnetic order after irradiation. 



Makarova. Han. Fsquinazi. da Silva. Kopelevich. Zakharova. and Sundqvistl (j200,^ 



studied the magnetism in photopolymerized fullerene films by MFM. The mag- 
netic force gradient measurements revealed that laser- and electron-beam ir- 
radiation of fullerene films produces magnetic images highly correlated to the 
topography. This correlation is expected since the polymerization shrinks the 
film material at the surface within a certain penetration depth, which depends 
on the radiation characteristics, producing topographic clusters with magnetic 
order. An example of the AFM and MFM images in a region illuminated in 
air with an energy of 2.6 eV and intensity of 200 mW/cm 2 is shown in Fig. 13. 
The high correlation between of magnetic grain formation and topography is 
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Fig. 13. Left side sketch shows the fullerene film with the two irradi- 
ated regions: the upper ellipse corresponds to the laser illuminated re- 
gion ( Makarova. Han. Esquinazi. da Silva. Kopelevich. Zakharova. and SundqvistL 



2003h and the lower rectangle to the region where 20 spots of 1.8 ^m diameter 



each where irradiated with a proton microbeam with a current of 500 pA and a 
fluence between 0.068 and 68 nC//mi 2 . In contrast to the spot irradiation in HOPG 
(see section 3.1.2) no change in the topography has been detected after irradiation 
of the spots in the fullerene film. All images correspond to an area of 5 x 5 |im 2 . 



well recognizable. In contrast to the laser illumination on Cgo film and pro- 
ton irradiation in HOPG, proton irradiation of spots of micrometer size on 
the fullerene film do not produce appreciable changes in the topography. The 
magnetic domains are however notable, see Fig. 13. 



These results obtained after irradiation of the thin films are unexpected and 
indicate either the hydrogen concentration in the sample or surrounding before 
irradiation is relevant and/or the defect concentration in the first micrometer 
from the sample surface and produced by the proton beam is not negligible. 
Because the thickness of the film is below 1 ixm most of the protons go through 
the material. 
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4 Annealing and Aging Effects 



We have measured with the SQUID some of the proton irradiated samples 
months after the irradiation. Depending on the sample, the enhancement of the 
magnetic order decreased or even vanished after more than 8 months leaving 
them at room temperature and in air. However, the magnetic signal remained 
for some of them even after one year. The reason for the different behavior is 
not known yet. In this section we shall summarize results that provide evidence 
for aging effects in the magnetic ordering, as well as annealing effects at high 
temperatures, measured by SQUID or MFM. Taking into ac count previous re- 
ports on the behavior and diffusion of hydrogen in graphite (|Atsumil . Eoolh as 



well as aging effects at room temperature in the magnetization of the fullerene 



Cfinli ^ ([Antonov. Bashkin. Khasanov. Moravskv. Morozov. Shulga. Ossipvan. and Ponvatovskv 



2002) one may expect to observe some time dependence in the magnetic re- 



sponse at the irradiated surface by MFM (or in the bulk magnetic moment 
by SQUID) if H is involved in the induced magnetism and a diffusion takes 
place. Two different measurements were done to test the possible influence of 
H diffusion on the magnetic signals. We measured different irradiated spots on 
HOPG and irradiated amorphous carbon film just after irradiation and after 
leaving the samples several months at room temperature. Other samples we 
annealed at high temperatures in vacuum or He atmosphere. 



Figure 14 shows the magnetic moment of a disordered carbon sample before, 
just after proton irradiation and after one year at room temperature (after 
irradiation). It is clearly seen that the ferromagnetic-like hysteresis vanishes 
after one year. 

Figures 15(a) and (b) shows the topography (top) and MFM (bottom) images 
of two spots and their surroundings created with 0.2 nC//xm 2 and 24.5 nC/ytmi 2 
on sample 1. The measurement was done one day after proton irradiation. 
Line scans of the topography (top) and MFM (bottom) images are shown in 
Fig. 15(c). As shown in this figure a clear enhancement in both topography and 
MFM signals at the irradiated region is clearly found. The enhancement of the 
topographic swelling increases with fluence with our irradiation conditions. For 
example, the peak height is ~ 8 nm (~ 225 nm) for a fluence of 0.2 nC//xm 2 
(24.5 nC//iin 2 ). However, there is no clear correlation between the increase 
in swelling height measured by AFM and the maximum phase shift change 
of the MFM signal at the spot within the fluence range used. As shown in 
Fig. 15(c) bottom, a magnetic structure within the spot produced with 24.5 
nC//im 2 appears, an indication for the existence of magnetic domains within 
the area of the irradiated spot. 

Figure 16 shows the same measurements as in Fig. 15 for the same two spots, 
but obtained 8 months later. After this time we found clear changes in the to- 
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Fig. 14. Magnetic moment (the diamagnetic background was subtracted) as a func- 
tion of applied field for the disordered carbon film CH577a after preparation by 
PLD (before irradiation (□)), after irradiation (•) and one year later, leaving the 
sample at room temperature (A). 

pography and in the phase shift of the MFM signal at and in the surroundings 
of the original microspots, see Fig. 16. There are three main changes to be 
noted: (1) The sharply white spots in MFM images (see the bottom Figs. 15(a) 
and (b)) on the irradiated area changed to black spots (see bottom Figs. 16 
(a) and (b)) and the maximum value of phase shift changes from ~ +0.7° 
to ~ —0.2°. It means that the direction of magnetization and its magnitude 
at the spots changed. (2) There appear clear and localized MFM signals be- 
sides the original spot signal (black spots in Figs. 16(a) and (b)). We note 
also that these "mirror" spots for both fluences are not clearly related to the 
topographic swellings as can be seen by comparing the line scan signals of to- 
pography and MFM of both spots (see Fig. 16(c)). There is also a clear change 
in the topography at the irradiated spots indicating a lattice relaxation. (3) 
The MFM images reveal that after 8 months there is a change not only at 
the spot position and its near neighborhood but also at its surroundings. This 
change is clearly observable in the vicinity of the spot made with a fluence 
24.5 nC/^m 2 where a well discernible domain pattern appears, see bottom 
Fig. 16(b). 



Hydrogen desorption in graphite was and still is a matter of research. Depend- 
ing on the characteristics of the H-trapping this desorption may or may not 
be accompanied by a lattice change. Therefore, it is not straightforward to 
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Fig. 15. Topography (top) and MFM (bottom) images (20 x 20 ^m 2 ) of sample 1 
for two spots created at 0.2 nC///m 2 (a) and 24.5 nC//jm 2 (b) and their surround- 
ings measured one day after irradiation. The tip-to-sample distance was 50 nm. (c) 
Corresponding line scans (white lines in (a) and (b)) of topography (top) and phase 
shift (bottom) images. 



conclude that the aging effects we observed are due only to the H-diffusion 
without any structural relaxation. Because we got evidence for a structural 
relaxation, this can influence the effective diffusion of H in the carbon ma- 
trix. Therefore, the effective activation energies for H diffusion in our system 
are not necessarily the same as, for example, those obtained by experimen- 
tal methods, usually at high tempe ratures, to study kinetics of diffusion of 
hydrogen in graphite (Atsumi, 20021 ). 



Certainly, one does not necessarily wait one year to see effects that may be 
related to hydrogen diffusion. It is quicker to anneal an irradiated sample 
at high enough temperatures and check for its influence. Figure 17 shows the 
MFM images obtained for the crosses corresponding to those of Fig. 10(a) and 
(c) measured at identical conditions after annealing the sample 2 hs at 1000 C 
in vacuum. Whereas no signal is detected for the cross in (b) ((c) in Fig. 10) 
magnetic domains are still identified after annealing in the surroundings of 
the cross (a). Compare this observation to that obtained after aging the spots 
at room temperature for several months shown in Fig. 16. 
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(a) (b) 

Fig. 17. The figures (a) and (b) show magnetic images of the same crosses shown in 
Fig. 10(a) and (c) measured under similar conditions with similar MFM tips after 
annealing them in vacuum at 1000 C for 2 hs. Whereas the crosses are not anymore 
visible in MFM, domain structure is still detected in (a). 
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5 Conclusion and Open Issues 



In this chapter we have reviewed the main observations after proton irradi- 
ation of carbon structures. The main results indicate that proton irradiation 
can trigger magnetic ordering in graphite. The results after proton irradiation 
leave no doubt that magnetic ordering exists in a carbon structure without 
the influence of magnetic ions. Neither the total amount of magnetic impuri- 
ties is sufficient to account for the measured magnetization nor the creation of 
magnetic spots in the micrometer range with the proton micro-beam can be 
understood based on magnetic metal-ion impurity concentration below 1 ppm 
as the PIXE results indicate. The results obtained up today do indicate that 
different parameters and sample states may play a role in determining the 
strength of the induced magnetic ordering. Broad irradiation (i.e. with a pro- 
ton beam of large diameter) appears to be less effective. The reason for this 
behavior will be studied in the future, but we speculate that it may be due 
to the produced defect density. Further experimental characterization (using 
the broad spectrum of methods in magnetism research) and sample prepara- 
tion studies are necessary to understand and stabilize the magnetic ordering 
found in carbon structures. The following issues should be clarified in the near 
future: 

(1) The role of H- atoms, implanted by irradiation as well as those already in 
the sample. 

(2) The contribution to the magnetic order from lattice defects produced by 
irradiation and their possible influence as H-trapping centers. 

(3) Dependence of the magnetic order on the induced type of defect. 

(4) The maximum achievable saturation magnetization in carbon structures. 

(5) The range of Curie temperature. 

(6) Influence of the ion current, fluence and energy of the irradiated particle 
on the magnetic order. 

(7) Influence of sample temperature. 

(8) Influence of the irradiation angle with respect to the crystallographic 
c-axis of graphite. 

(9) Influence of ion irradiation (other than proton) on the magnetism of 
carbon structures. 

(10) The effective magnetic moment of magnetic impurities in graphite as well 
as in disordered carbon structures. 

An answer to part of these open questions will take several years of research. 

Summarizing briefly recently published theoretical work, which is correlated 
to the influence of hydrogen in carbon magnetism we note that: 

(1) Hydrogenated nanographite can have a spontaneous magnetization due to 
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(2) 



(3) 



(4) 



(5) 



differ ent numbers of mono- and dihydrogenated carbon atoms ((Kusakabe and Maruvama . 
20031 ). Theoretical work using local-spin-density approximation calcu- 



lates the spin polarization of the graphite bands when different atoms 
as hydrogen, fluor or oxygen are added on a graphene layer and predicts 
magnetic orderin g with fully or partially sp i n-pola rized flat band, upon 
the added atom ( Maruvama and Kusakaba 2004 ). Theoretical simula- 
tions indicate that a magnetic ordered state should be easier to achieve 
with a mixture of carbon-hydrogen bonds than fluorinated nanograph ite 



( Maruvama. Kusakabe. Tsunevuki. Akagi. Yoshimoto. and Yamauchil . 12004). 
Although the calculations were done attaching the hydrogen at the car- 
bon atoms at the edges of a graphene layer, one may speculate that 
hydrogen can trigger the sp 2 -sp 3 unbalance promoting magnetic ordering 
in different carbon structures. 

Recently published work indicates that hydrogen modifies substantially 
the electronic structure of graphite around it. According to STM/AFM 
measurements ( Ruffieu x. Groning. Schwaller. Schlapbach. and Groning , 2000f ) 
a single H-atom interacting with a graphite surface modifies the electronic 
structure over a distance of 20 to 25 lattice constants. 



Muon spin rotation/relaxation experiments ([Chakhalian. Kiefl. Dunsiger. MacFarlane. Miller. Soe 

2002) indicate that a positive muon in graphite trigger s a local magnetic 

mom ent around it. Recently published theoretical work ((Duplock. Scheffler. and Lindan , 



20041 ) supports this conclusion. It also indicates that upon the type of de- 



fect in a graphene layer, hydrogen may not trigger any magnetic ordering. 

The magnetic moment and diffusion of adatom defects in a graphite sheet 

were s tudied bv lLehtinen. Foster. Avuela. Krasheninnikov. Nordlund. and Nieminen 

(2003) . The results of a full spin-polarized density functional theory in- 
dicate that these defects may have a magnetic moment of about 0.5 ur. 
The ca lculations of lLehtinen. Foster. Ma. Krasheninnikov. and Nieminen 

(2004) indicate that if an hydrogen encounters an empty vacancy, then it 
compensates the dangling bond and the magnetic moment of the vacancy 
vanishes. However, if a vacancy is saturated by an hydrogen atom, a sec- 
ond hydrogen atom will bond to the other side of the vacancy having a 
magnetic moment of 1.2 \xb localized o n the dangling sp 2 bond. 



New theoretical and experimental work ([Chan. Montanari. Gale. Bennington. Tavlor. and Harriso 
12004 ) suggests that hydrogen may play also an important role on the 
magnetic ordering found in fullerenes. 
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